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IMPROVED PREDICTION OF LAMINAR LEADING EDGE SEPARATION 
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McDonnell Aircraft Con^any 
McDonnell Douglas Corporation 


I SUMMARY 

I This research was conducted In order to provide a definite criterion for 

i the prediction of the bubble burst on airfoils typical of those used for fight- 

I er wings. The approach taken was to correlate existing airfoil bubble burst 

I data using various parameters at the laminar separation point. The method 

I due to Weber was modified to provide a continuous analytic solution 

I for the velocity distribution around the airfoil leading edge. This 

proved to be as accurate as any available. Coupling the modified Weber 
method with the Stratford laminar separation prediction method leads to a 
I universal chart giving the conditions at separation as a function of stagna- 

tion location and leading edge radius. Application of the cond>lned method to 
available two-dimensional airfoil data resulted In an empirical criterion 
presenting *:he limiting local velocity gradient at separation as a function of 
the boundary layer momentum thickness at separation for bubble burst. The 
correlation leads as well to the qualitative explanation of two types of lam- 
inar stall: thin airfoil and leading edge. The validity of the correlation 

Is demonstrated by predicting the lift coefficient and angle of attack for 
stall on airfoils with leading edge or trailing edge flaps. 


INTRODUCTION 

Airfoil laminar bubble burst Is the cause of excessive drag due to lift 
on very thin airfoils and of maximum lift stalling on moderately thin airfoils 
It has been the subject of considerable experimental and theoretical research. 
However, efforts to apply previous research to the subject of maneuvering drag 
due to lift on swept thin wings typical of transonic fighter aircraft proved 
unrewarding. The present research was undertaken to provide the basis of a 
technique to predict laminar bubble burst on such w.'ngs. 

As pointed out by Goradla and Lyman (Reference 1) bubble burst Is quite 
likely to be a function of conditions at the laminar separation point. The 
problem therefore consists of determining those conditions and correlating 
the experimental data. 


NOMENCLATURE 


A 

A(x) 
B 


Average value of A(x) 

Velocity ratio due to thickness at a ■ 0 
Average value of B(x) 
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®(*) Velocity ratio due to thickness at angle of attack 

Lift coefficient 

g One half the leading edge radius divided by chord 

Rg Freestream Reynolds number based on chord 

^0 Local Reynolds number based on momentum thickness 

s Arc length along the airfoil surface, divided by chord 

u The ratio of local velocity to freestream velocity 

Up Peak value of the velocity ratio 

Ug Separation value of the velocity ratio 

V Local velocity 

Freestream velocity 

X Distance from the nose along the chord divided by chord 

x^ Stagnation location 

a Angle of attack 

LAMINAR SEPARATION 

As discussed in Reference 2, the theory of Stratford, Reference 3, is 
very accurate. But, like most other methods it requires accurate, smooth 
longitudinal pressure gradients. Utilization of experimental pressures with 
attendant imprecision of measured pressures and pressure tap location proved 
impossible. Using the currently available powerful finite difference theories 
presented equal difficulties due to the inherent discontinuous pressure gra- 
dient. 

A new method was derived based on that of Weber, Reference 4, While not 
mathematically rigorous, this method has proved to be of very high precision 
near the airfoil nose both in pressure level and in pressure gradient. It 
has been derived for the general airfoil shape, but is best illustrated for the 
uncambered airfoil. The local velocity is given by 


V 

00 


U 


A(x) Cos a + B(x) Sin a 



ds 

dx 


The functions A(x) and B(x) are functi''ns of the thickness distribution, but 
for most airfoils are nearly independent of x. Furthermore, the surface slope 
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function can be approximated by a parabola. Thus, the velocity aai gradient 
are readily found by 


u 


A 


cos B Sin 

V^+ X 


a 


du 

ds 


i( A g Cos g - (+B Sin a) ] 


(x + g)‘ 


The stagnation (V = 0) location Is therefore. 


X 2 2 

s tan g 

g ” ^ g 

It Is convenient to give the local velocity in terms of the stagnation 
location, rather than angle of attack. It Is apparent that the airfoil velo- 
city distributions become universal functions of x/g and Xg/g as shown in 
Figure 1. It is also apparent that certain key events are also functions of 
these parameters as shown. In particular, the laminar separation point and 
conditions at that point are direct functions of Xg/g. 

The Stratford method can be applied to the non-dimensional velocity dis- 
tributions to derive a universal chart for the velocity, position, and gra- 
dient at separation. This is given in Figure 2. The figure can be used for 
any airfoil once the stagnation point is known. This formulation is particu- 
larly useful in that the correct non-dimensionalizing paramecerr become obvious. 

BUBBLE BURST CORRELATION 

The available two-dimensional airfoil data was analyzed using this method 
to determine boundary layer and velocity gradient parameters at the laminar 
separation point at or near the angle of attack for laminar bubble burst. The 
best correlation obtained i? shown in Figure 3. The data scatter results from 
a + 1/2 degree angle of attack precision in the bubble burst angle of attack 
data. It should be noted that velocity gradient parameter ia a 3 /7th power 
function of the Reynolds number parameter in keeping with the general postu- 
lation that bubble burst is a failure of the turbulence level to re-atcach. 

Figure 3 ia a key element of airfoil analysis but is not the only factor. 

If the Reynolds number is high enough, transitions will occur prior to laminar 
separation, precluding bubble formation. There is no wind tunnel data avail- 
able at high enough Reynolds number to demonstrate such a phenomenon. But, 
on the other hand, there is considerable data at very low Reynolds numbers to 
indicate what has been called long bubble stall. Presumably this occurs when 
the separated laminar boundary layer does not immediately transition to turbulent. 
It is therefore unable to re-attach and the so-called short bubble does not 

form. It is generally thought that this occurs for R* < 125. 
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Post stall analyses are very difficult because of the total reliance on un- 
reliable experimental pressure distributions. However » there is good evidence 
to lead to the conclusion that as the short laminar bubble bursts » the pres- 
sure distribution must adjust, with Increasing angle of attack, so that condi- 
tions for laminar bubble burst are maintained. In other words, the bubble 
burse boundary applies not only for the burst angle of attack, but for post 
burst angles as well. 

Figure 4 shows the sequence of events and the relation to two well kno%m 
phenomena: leading edge stall, and thin airfoil stall. For instance, airfoil 
A is typical of moderately thin 'irfolls (i.e., 64A010). At very low angles 
attack, laminar separation does not occur near the leading edge, but aft of 
maximum thickness (if at all). As angle of attack increases the separation 
moves to the airfoil nose and a short bubble forms. As angle of attack is 
increased up to the burst boundary, a typical short bubble is obtained with 
kittle impact on the overall forces and moments. Further increases in an,'$le 
of attack cause a drastic reduction in lift as the laminar separation condition 
moves down the boundary* 

In the case of airfoil B, typical of thinner airfoils (i.e., 64A006), the 
same sequence of events occurs. However as the bubble bursts, and as angle 
*^ttack is increased, conditions at laminar separation move below R0 = 125. 

The separated layer does not transist and the long bubble is formed. Lift con- 
tinues to increase, but at a reduced rate. 

MAXIMUM LIFT PREDICTION 

While this research was aimed primarily at drag prediction at high angles 
of attack, it does form the basis for prediction of maximum lift for airfoils 
of the "leading edge stall" type. This is xllustrated in Figure 5. The ef- 
fects of a leading edge flap are shown on a 64A010 airfoil, and of a trailing 
edge flap on a NACA 0009 airfoil. In the former case, it was noted during the 
test that separation occurred at the hingeline for flap angles in excess of 
15 . Up to this angle, the current method predicts both Cl^, and angle of 
attack for CL^ax very well. 

In the case of the trailing edge flap, no attempt was made to adjust the 
predicted airfoil nose pressure distribution for viscous effects on the flap 
itself. Even so, the angle for CL„,ax predicted and the prediction for 

^Lmax consistent with the approximation. 

CONCLUDING REMARKS 

Rational approximation to airfoil theoretical pressure has yielded an 
accurate method for prediction of conditions at laminar separation and given 
insight into a universal correlation of the bubble burst phenomenon. This, 

in turn, has led to a basic understanding of thin airfoil and leading edge 
s tisH • 
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Figure 1.- T 3 rplcal nondlmensional velocity profiles. 



’‘stages 


Figure 2.- Conditions at the laminar separation. 
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Figure 3.- 
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Correlation of parameters affecting laminar bubble burst 



Figure 4.- Leading-edge stall development. 










